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Investigation  of  Closed  Loop  Adaptive  Optics  with  the 
Deformable  Mirror  not  in  Pupil  -  Part  2:  Theory 
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ABSTRACT 

This  is  the  second  of  two  papers  discussing  aspects  of  placing  the  deformable  mirror  in  a  location  not  conjugate 
to  the  pupil  plane  of  a  telescope.  The  Starfire  Optical  Range,  Air  Force  Research  Laboratory  is  in  the  process 
of  developing  a  high  efficiency  adaptive  optic  (AO)  system  for  its  3.5m  telescope.  The  objective  is  to  achieve 
maximum  optical  throughput  while  maintaining  diffraction  limited  performance.  This  can  be  acheived  by  mini¬ 
mizing  the  number  of  optical  elements  which  then  implies  placing  the  AO  system  on  gimbal.  Doing  so  forces  the 
deformable  mirror  (DM)  to  not  be  conjugate  to  the  pupil.  However,  placing  the  DM  so  that  it  is  not  conjugate 
to  the  pupil  results  in  a  degradation  in  Strehl.  Wave  optics  simulation  has  shown  a  7-15%  degradation  in  Strehl. 
In  the  previous  paper,  we  experimentally  measured  this  degradation.  In  this  paper,  we  determine  the  cause  of 
the  degradation.  Then  we  determine  a  means  to  mitigate  this  degradation  in  a  manner  that  is  implementable  on 
conventional  AO  sytems.  Simulation  results  are  presented  to  validate  the  veracity  of  the  mitigation  technique. 

Keywords:  Adaptive  Optics,  High  Efficiency  AO 

1.  INTRODUCTION 

Many  astronomical  observatories  are  designing  their  AO  systems  with  the  deformable  mirror  not  conjugate  to  the 
pupil.  This  is  done  for  a  variety  of  different  reasons  ranging  from  the  need  to  perform  multi-conjugate  adaptive 
optics  (MCAO)  to  the  need  to  minimize  the  number  of  optical  surfaces  in  order  to  improve  throughput.  An 
innovative  approach  used  by  some  telescopes  converts  the  secondary  mirror  into  a  deformable  mirror  (DM).  At 
wavelengths  longer  than  2fj,m  where  there  is  a  need  to  limit  the  number  of  warm  surfaces,  an  optical  design 
with  an  adaptive  secondary  minimizes  surfaces  thereby  maximizing  sensitivity  making  it  particularly  suitable  for 
thermal  infrared  astronomy.  For  infrared  astronomy  an  adaptive  secondary  is  also  attractive  for  MCAO  since 
the  secondary  is  roughly  conjugate  to  the  strong  ground  layer. ^ A 

The  6.5m  Multiple  Mirror  Telescope  (MMT)  at  Mt.  Hopkins  in  Southern  Arizona  is  the  first  in  the  world  to 
use  an  adaptive  secondary^A  fQj-  wavefront  sensor  correction.  Placing  the  DM  at  secondary  in  this  case  allows 
for  higher  throughput,  lower  emissivity,  and  a  simpler  optical  setup.  Miller  et  al."^  report  the  measured  Strehl  to 
be  40  %  of  the  diffraction  limited  peak. 

The  Large  Binocular  Telescope  (LBT)  on  Mt.  Graham,  Arizona  has  two  8.4m  primary  mirrors  and  two 
0.91m  Gregorian  secondary  mirrors.  The  secondary  mirrors  are  the  deformable  mirrors  for  the  LBT  adaptive 
optics  system.^  The  DMs  for  the  LBT  were  designed  based  on  the  successful  MMT  adaptive  secondary  design. 
The  Very  Large  Telescope  (VLT)  at  Paranal,  Chile  will  upgrade  one  of  the  8m  telescopes  to  include  a  deformable 
secondary  mirror  (DSM).® 

The  Starfire  Optical  Range  (SOR)  is  also  developing  a  high-efficiency  AO  system  for  its  3.5m  telescope.  One 
proposed  configuration  replaces  the  turning  flat  M4  with  the  DM.  In  this  geometry,  the  DM  is  lAlm  (in  output 
space)  behind  the  3.5m  primary  mirror.  Waveoptic  simulation  of  this  geometry  by  Foucault®  indicated  that  this 
geometry  leads  to  a  7-15%  degradation  in  Strehl. 
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This  research  seeks  to  replicate  this  degradation  experimentally,  i.e  verify  that  the  degradation  is  not  a 
simulation  anomaly,  to  find  the  root  cause  of  the  degradation,  and  then  to  determine  if  a  means  to  mitigate 
the  degradation  was  possible.  As  such,  our  research  consisted  of  two  parts,  experimental  and  theoretical.  Ex¬ 
perimentally,  we  built  an  AO  testbed  of  the  proposed  DM-not-in-pupil  (DNiP)  geometry^  and  conducted  a 
set  of  experimental  tests.  This  experiment  verified  the  existance  of  the  degradation  and  we  measured  a  7-9% 
degradation  in  Strehl.  These  results  are  presented  in  Part  1  of  this  two  part  series. 

In  this  paper,  we  discuss  the  root  cause  of  the  degradation  and  also  present  a  means  to  mitigate  the  degra¬ 
dation  in  such  a  way  that  it  is  implementable  on  conventional  AO  systems.  To  this  end,  Section  2  details  SOR’s 
DM-Not-in-Pupil  (DNiP)  geometry  under  which  the  problem  arises.  Section  3  determines  the  root  cuase  of  the 
degradation.  Section  4  calculates  a  closed  form  algebraic  solution.  Section  5  proposes  a  theoretical  description 
for  the  DNiP  degradation.  Section  6  presents  a  mitigation  techniuqe,  and  Section  7  concludes  our  findings  on 
closed  loop  Strehl  performance  with  DNiP. 


2.  BACKGROUND 

A  conventional  AO  system  places  the  DM  conjugate  to  the  pupil.  Doing  so  causes  rays  at  the  DM,  aside  from 
magnification,  to  have  the  same  height  and  angle  as  rays  at  the  pupil.  In  this  configuration,  the  science  and 
beacon  beams  at  both  the  pupil  and  DM  have  the  same  geometry.  So  when  the  DM  uses  the  beacon  beam 
measured  at  the  DM-plane  to  correct  the  science  beam  in  the  DM-plane,  the  AO  correction  is  always  proper. 
Placing  the  DM  so  that  it  is  not  conjugate  to  the  pupil  does  not  necessarily  cause  the  correction  to  be  in  error;  for 
instance  when  natural  guidestar  AO  is  used  and  the  pupil-DM  separation’s  Fresnel  number  is  high,  the  correction 
will  be  proper. 

However,  the  geometry  shown  in  Figure  1  does  induce  a  degradation  in  Strehl.  In  this  geometry,  (I)  the 


Figure  1.  DNiP  geometry.  Left:  Mf  has  been  replaced  by  a  DM  causing  a  displacement  of  approximately  IKm 
in  outputspace.  Right:  the  science  object  is  at  infinity  and  the  sodium  guidestar  at  approximately  90Km.  These  two 
conjugates  lead  to  the  sodium  beam  diverging  after  the  telescope  removes  the  focus  of  the  science  object. 


telescope  is  focused  on  the  science  object,  so  the  science  beam  is  (stochastically)  planar  after  the  primary¬ 
secondary  optical  elements,  (2)  the  DM  is  placed  non-conjugate  to  the  pupil,  and  (3)  the  sodium  beacon  is 
generated  at  approximately  90Km  causing  the  sodium  beacon  to  have  a  residual  (de)focus  after  the  primary¬ 
secondary  optical  elements.  This  is  the  geometry  of  one  proposed  design  for  the  SOR  high  efficiency  AO  system 
which  replaces  M4  by  a  DM  causing  a  displacement  of  approximately  IKm  in  outputspace.  This  is  the  design 
for  which  a  7-15%  degradation  is  Strehl  was  noted  and  the  design  to  be  studied  here. 
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3.  DETERMINING  THE  ROOT  CAUSE  OF  THE  DEGRADATION 


There  are  well  known  causes  for  degradation  in  Strehl  and  this  degradation  was  initially  thought  by  the  commu¬ 
nity  to  be  caused  by  focus  anisoplanatism.® 


3.1.  Focus  Anisoplanatism 

Focus  anisoplanatism  arises  with  LGS  AO  because  the  LGS  is  at  a  finite  altitude  with  respect  to  the  science  object 
(see  the  left  plot  in  Figure  1).  Because  the  beacon  is  at  a  finite  altitude,  it  does  not  sample  the  full  atmosphere 
experienced  by  the  science  beam.  This  improperly  sampled  atmosphere  is  then  used  as  the  correction  for  the  full 
atmosphere  causing  a  residual  phase  variance  on  the  science  beam  even  if  the  measured  phase  is  fully  corrected. 

This  effect  is  well  known  and  has  been  thoroughly  studied.  Following  Sasiela,®  the  phase  variance  caused  by 
focus  anisoplanatism  for  a  spherical  wave  is  given  by 

D5/3 

cr2=o.5fc^—  M5/3_, 

with  ^5/3=/  dz 

Jo 

where  D  is  the  pupil  diameter,  ko  is  the  wave  number,  L  is  the  distance  to  the  guide  star,  and  is  the 
atmospheric  structure  function,  /X5/3_  is  the  (5/3)*^  atmospheric  moment,  the  subscript  _  denotes  integration 
from  the  guidestar  down,  and  the  science  signal  is  at  infinity.  Throughout  this  analysis  is  given  by  HV-21 
(commonly  known  as  RV5/7),  i.e. 


(1) 

(2) 


C'2  =  0.00594^  (10"®/i)^°eTOw  +  2.710"^®eT5ro  +  Ae TOT  (3) 

where  h  is  the  altitude  in  meters,  is  in  A  =  1.710“^"^,  and  W  =  21  a  pseudo-wind  term. 

Note  that  if  the  disturbance  is  solely  in  the  pupil  then  C^{h)  =  C^{0)6{h)  which  implies  that  =  0,  which 
implies  S'  =  1,  i.e.  there  is  no  focus  anisoplanatism  regardless  of  the  value  of  AL.  Hence,  conventional  wisdom 
is  incorrect. 

But  perhaps  subtleties  occur  when  is  distributed.  So,  assuming  a  distributed  atmosphere  Equation  1  can 
be  rewritten  as 


cr^  =  0.5U 


D 

L  +  AL 


5/3 


L  +  AL  —  z 
L  +  AL 


:)^/^dz 


where  AL  is  the  distance  out  of  the  pupil.  Using  the  standard  approximation,  the  degradation  in  Strehl  given 
by  this  phase  variance  is 

S'  ~  e-*"'  (4) 

and  so  the  degradation  caused  by  putting  the  DNiP  geometry  is  given  by  Hence,  the  change  in  focus 

anisoplanatism  due  to  the  DNiP  geometry  causes  a  decrease  in  Strehl  of  0.49%  when  A  L  is  increased  to  1  km. 
This  is  plotted  in  Figure  2.  Hence,  there  is  less  than  1%  additional  degradation  due  to  focus  anisoplanatism 
when  the  DM  is  placed  in  the  DNiP  position  for  a  realistic  atmosphere.  This  conclusion  is  unchanged  if  the 
piston  and  tilt  focus  anisoplanatic  phase  variance  is  used.  So,  the  degradation  seen  in  simulation  and  experiment 
are  not  caused  by  focus  anisoplanatism. 


3.2.  ro,  6q,  cj^,  and  Greenwood  frequency 

There  are  a  standard  list  of  stochastic  parameters  that  cause  degradation  in  Strehl.  It  was  thought  that  perhaps 
the  degradation  in  Strehl  was  caused  by  changes  in  these  parameters  and  this  caused  the  degradation.  So, 
studying  how  these  parameters  change  with  changes  in  AL  is  necessary.  Specifically,  this  subsection  modifies  the 
expressions  for  ro,  Rytov,  and  Oq  to  account  for  the  DM-not-in-the-pupil  (DNiP)  case.  From  physical  grounds, 
the  DNiP  degradation  occurs  even  when  the  science  object  is  a  star,  so  while  the  Greenwood  frequency  causes 
AO  system  degradation,  it  plays  no  part  in  the  DNiP  degradation  and  will  not  be  discussed  further. 
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Strehl 


AL  [km] 


Figure  2.  Degradation  in  Strehl  caused  by  focus  anisoplanatism  as  a  function  of  AL 


3.2.1.  Ftied’s  Parameter,  Rytov,  and  the  Isoplanatic  Patch 

Again  following  Sasiela 


ro  =  (0.423A:^^o)”^/^  =  |^0.423A:g  ^  Cl{z)  dzj 

pL 

al  =  0.563A:o"/V5/6  =  0.291A:o"/®  /  Cl{z){^zf/^  dz 

6»o  =  2.91fcoM5/3  =  2.91A:o  f  Cl{z){^zfl‘^  dz 


(5) 

(6) 
(7) 


with  7  =  for  a  laser  guidestar  geometry. 

To  modify  the  equations  for  the  DNiP  geometry  with  a  spherical  beam,  L  is  replaced  with  L  +  AL  and  the 
integration  is  from  —  AL  to  L.  So  for  a  spherical  beam  propagating  from  L  to  0,  7  =  ■  Then,  the 

expressions  become: 


ro=^0.423fco^  Cl{z)d^ 


-3/5 


cr^  =  0.291fco 


7/6 


6»o  =  2.91A:, 


' 

J-AL  L-\-AL 

f  /-i2i  x/L-\-AL —  Z  x5/3  7 

r  +  AL 


But  C^iz)  =  0  when  z  <=  0,  so 


ro  =  ^0.423A:o  J  C^(z)  dz'^ 

=  0.291fco"/®  c2(z)(^^±^4^z)5/6  dz 

Jo 


L  +  AL 
,  L  +  AL  —  z 


00  =  2.91fc2  ^ 


(8) 

(9) 

(10) 


(11) 

(12) 

(13) 


Equation  11  demonstrates  that  rg  does  not  change  with  AL.  So,  changes  in  rg  cannot  be  responsible  for  the 
DNiP  degradation. 
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Equations  12  and  13  have  a  form  such  that  AL  only  appears  when  added  to  L.  So  for  AL  <  L,  changes  in 
rg  and  Rytov  should  be  small.  This  is  indeed  the  case;  note  as  AL  increases  from  0  to  5  km,  increases  by 
0.32%,  while  9o  decreases  by  0.59  %.  The  0.32%  increase  in  scintillation  corresponds  to  a  multiplicative  decrease 
in  Strehl  of 

The  decrease  in  the  isoplanatic  patch  causes  a  a  change  in  Strehl  of 


with  e  =  0.0059  when  AL  =  5Km  and  the  rightmost  relation  holding  for  any  6  <  Oq  and  AL  <  bKm.  So,  the 
DNiP  degradation  is  not  caused  by  changes  in  scintillation  or  the  isoplanatic  patch. 

These  results  are  plotted  in  Figures  3,  4,  and  5.  In  each  plot,  the  beacon  is  at  100  km,  A  =  and  the 

model  used  is  lx  the  HV-21  model. 

Section  Summary  This  section  has  demonstrated  that  the  cause  of  the  DNiP  degradation  is  not  one  of  the 
customary  degradation  effects.  By  process  of  elimination,  the  only  remaining  physically  varying  quantity  is  the 
variation  in  magnification  between  the  pupil  and  DNiP  planes.  We  propose  that  this  must  be  the  cause  of  the 
degradation.  This  result  was  reached  independently  by  Rhoadarmer.^^  In  keeping  with  standard  notation  in  the 
AO  field,  we  have  named  this  degradation  “magnification  anisoplanatism” .  The  next  few  sections  work  through 
this  hypothesis. 
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Figure  3.  ro  as  a  function  of  AL 
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Figure  4.  Scintillation  as  a  function  of  A  L 
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Figure  5.  9q  as  a  function  of  AL 

4.  AN  EXACT  EXPRESSION  VIA  FRESNEL  PROPAGATION 

In  order  to  gain  a  basic  understanding  of  the  problem,  this  section  studies  the  formalism  that  propagates  the 
light  from  the  pupil  plane  to  the  DNiP  plane.  These  propagations  are  well  known  to  be  Fresnel  propagations. 

4.1.  Fresnel  Propagation  from  the  Pupil  to  the  DNiP  Plane 

In  general,  a  field  at  the  (x,y)  plane  propagated  from  the  plane  is  described  as: 

ikz  foo  POO 

U{x,y)  =  —  dU  (15) 

J  —  oo  J  —  oo 

where  A  is  the  wavelength,  A:  =  ^  is  the  wavenumber,  and  U (^,  ry)  is  the  field  at  the  source  plane.  The  effect  of 
a  lens  is  included  by  multiplying  Equation  15  by 

b(a;,y)  =  (16) 

where  /  is  the  focal  length  of  the  lens.  In  the  DNiP  case,  the  focal  length  of  the  lens  at  the  pupil  /  =  Zgd-  So 
for  the  target,  the  focus  error  found  in  the  phase  at  '>  is  removed  with  the  lens  term  ^ 

and  the  resulting  equation  for  the  target  signal  immediately  after  the  lens  is 

pikZsci 

=  - -  (17) 

tAZsci 

Similarly,  the  equation  for  the  signal  from  the  beacon  immediately  after  the  lens  is  given  by 


Ui{^,v)  =  e 


iXzh 


but  in  this  case,  the  focus  term  is  not  completely  canceled  by  the  lens. 


So  the  residual  focus  term  is  given  by  e  ^sai' ,  The  problem  can  be  simplified  further  if  we  assume 

that  the  signal  from  the  target  and  the  signal  from  the  beacon  have  uniform  intensity  at  the  aperture,  with 
the  added  residual  focus  term  for  the  signal  from  the  beacon.  This  is  a  good  approximation  for  astronomical 
imaging,  where  the  Rytov  is  less  than  0.1.  Then,  Usci{^,'n)  =  1  inside  the  aperture  with  diameter  D  and 
=  e  ^  '  ''■‘b  ‘sai’  inside  the  aperture  with  diameter  D.  These  can  be  propagated  to  the  DM  plane 
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(x,y)  {zdm  from  the  pupil)  to  determine  the  effect  of  propagating  the  focus  error  on  the  phase  variance  at  the 
DM  plane. 


Usci{x,y) 


Ub{x,y) 


gikzdm 


d^drj  e 


pupil 


d^drj  e 


pupil 


(20) 

(21) 


These  two  equations  can  be  compared  for  various  scenarios  but  what  we  really  want  is  perfect  compensation 
of  the  beacon,  i.e.  calculate  the  phase  from  the  beacon  at  the  DM  plane  and  subtract  this  from  the  phase  from 
the  target  signal  at  the  DM  plane.  The  Strehl  degradation  is  determined  by  comparing  the  variance  of  the 
residual  phase  at  the  DM  plane  to  the  variance  of  the  residual  phase  at  the  pupil  plane.  To  get  this  expression, 
assume  the  WFS  is  focused  at  the  beacon.  Then,  the  beacon  beam  has  no  focus  error.  Add  the  propagated 
residual  focus  error  (making  the  beam  spherical  again,  with  a  source  at  the  target)  and  back  propagate  to  the 
pupil  plane  to  re-create  the  expanding  beam  at  the  pupil.  Then,  remove  the  residual  focus  error  to  collimate 
the  beam  and  forward  propagate  to  calculate  what  a  collimated  field  from  the  beacon  would  look  like  in  the 
DM/WFS  plane.  This  is  done  in  Appendix  A  yielding 


ik 


Arg{U{x',y'))  =  Arg  —  (x'^-ky'^) 

■y' 


WFS 


dxdyUmeas{x,y)e 


ik{{y/ax-\-  - 


^  +  (V“!/+ 


(22) 


Section  Summary  With  this  expression,  we  have  shown  that  it  is  possible  to  determine  the  correction  needed 
in  the  pupil  plane  based  on  the  measured  data  -  -  the  magnification  hypothesis  appears  as  scaling  of  the  coor¬ 
dinates  on  the  right  hand  side  of  the  equation.  However,  while  this  is  a  closed  form  solution,  it  is  not  easily 
implementable  on  realtime  AO  systems. 

5.  PHYSICAL  DESCRIPTION  OF  THE  DEGRADATION 

Assuming  magnification  is  the  cause  of  the  degradation,  it  is  possible  to  construct  an  expression  which  corrects 
for  the  degradation  with  an  additive  term;  the  additive  term  allows  the  correction  to  be  implemented  within 
existing  AO  control  paradigms.  That  is,  by  using  the  gradient  reference  to  mitigate  the  degradation,  the  DM 
and  WFS  remain  both  conjugate.  In  this  manner,  the  DM- WFS  pair  remain  a  null  seeking  device  with  the  null 
being  a  full  correction  of  the  science  beam,  not  the  beacon  beam. 

Notationally,  (j){r)  is  used  throughout  as  the  science  object’s  phase  in  the  pupil,  (j){ar)  as  the  beacon  phase 
measured  in  the  DM- WFS  plane,  and  a  as  the  magnification,  a  has  been  calculated  to  be  «  1.06  for  the  proposed 
SOR  high  efficiency  AO  system,  but  for  ease  of  calculation  and  in  order  to  place  an  upperbound  on  performance 
a  =  1.1  is  used  throughout.  Also,  for  simplicity,  the  science  and  beacon  wavelengths  are  assumed  to  be  the  same. 

The  degradation  in  Strehl  is  approximated  by 

•SoNiP  =  (23) 

with  cr^  the  variance  in  phase  caused  by  a  non-unity  magnification  between  the  science  and  beacon  beams.  This 
phase  variance  is  given  by 

^ (jy  — ^  (^science  ^beacon)  ^  ^  ^science  ^beacon  ^ 

and  assuming  non-unity  magnification 
=<  ((/)(r)  —  (j){ar))‘^  >  —  <  (l>{r)  —  (j){ar) 

=<  Wr)  -  f  >  -(<  >f  .^4) 

and  because  the  atmosphere  is  a  zero  mean  process  ' 

=<  {4>{r)  —  (j){ar)Y  > 

=<  (A0(r))2  > 

where  A^(r)  =  (j){r)  —  4>{ar) 
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Then  assuming  a  coarse  grain  derivative 


<  (</,'(r))2  > 

where  5r  =  ar  —  r  \s  assumed  to  be  small  and  (j)'  {r)  is  the  derivative  of  the  phase  at  r. 
The  degradation  is  in  Strehl  due  to  the  DNiP  geometry  is  approximated  by 


•^DNiP  =  e 


with  (T^/  the  variance  in  the  derivative  of  the  phase. 


(25) 


(26) 


6.  MITIGATION  OF  THE  DEGRADATION 

Beginning  with  the  results  of  the  previous  section,  i.e.  assuming  magnification  is  the  cause  of  the  degradation  it 
is  trivial  to  write  a  formal  solution  that  mitigates  the  degradation,  i.e. 

(j){r)  =  (j)[ar)  —  A(j){r)  (27) 


with  <j>{r)  the  science  object’s  phase  in  the  pupil,  (j>{ar)  the  beacon  phase  in  the  DM- WPS  plane,  and  A(/)(r)  an 
additive  correction.  Then,  if  A(j>{r)  is  added  to  the  gradient  reference,  the  source  phase  is  written  on  the  DM, 
the  degradation  is  corrected,  and  formally  the  solution  is  in  hand. 

But  in  the  real  world,  A(j){r)  is  unknown.  So  considering  astronomical  imaging,  i.e.  <  0.2,  the  phase  is 
continuous.  As  such  (()(r)  can  be  expanded  in  a  Taylor’s  series  expansion  about  ar,^^  i.e. 


0(r)  =  ^<(.«| 


and 


Rearranging  terms  leads  to 


(j){r)  =  (f){ar)  +  {ar){ar  -  r)  +  X! 


k=2 


A(j){r)  =  (j){ar)  —  4>{r) 

=  —(j)^^\ar)  {a  —  l)r  —  S{r) 

with  £{r)  =  1^  error  term.  Then  substituting  into  Equation  27 


(28) 


(29) 


(30) 


(j){r)  =  (j)[ar)  +  <j)^^\ar)  (a  —  l)r  +  S{r) 
and  dropping  the  error  term,  we  obtain  this  section’s  main  result 

(j){r)  «  <p{ar)  +  (j)^^\ar)  {a  —  l)r 


(31) 


(32) 


Equation  32  is  the  means  to  correct  the  beacon  beam  phase  in  the  DNiP  plane  to  give  a  correct  AO  correction 
to  the  science  beam  in  the  pupil  plane. 


Estimation  of  the  Error  The  error  associated  with  this  correction  can  be  obtained  via  5(r)  =  |r 

Since  (j)  is  continuus  and  smooth  (recall  cr^  <  0.2),  6^  — >  0  as  fc  ^  oo  with  bk  the  term  in  £{r),  i.e. 

bk  =  .  This  allows  us  to  keep  a  finite  number  of  terms  with  confidence.  In  the  physical  system, 

(()(r)’s  maximum  rate  of  fluctuation  will  be  governed  by  Fried’s  parameter.  In  the  measured  and  applied  phase, 
the  maximum  rate  of  fluctuation  will  be  governed  by  the  DM’s  actuator  spacing.  In  either  case,  the  fluctuation 
fixes  the  maximum  order  of  a  polynomial  required  to  fit  the  phase.  This  sets  a  bound  on  k  above  which  the 


error  is  negligable  for  any  order  of  derivative.  For  instance,  consider  an  abnormally  high  magnification,  a  =  1.1, 
then  partially  determines  the  rate  of  convergence  for  small  k  and  fixes  it  for  large  k  yielding 

for  various  radii  the  following  errors 


r  =  0 

r  (0.1r)'‘  1 

1  fc'  J 

^  fc=2,“ 

=  {0,0,0,...} 

r  =  f  =  0.875 

r  (0.1r)'‘  1 

L  J 

fc=2,“ 

=  (3.8  X  10-3, 1.1  X  10-4, 2.4  X  10-®, . 

(33) 

r  =  f  =  1.75 

-1 

r  (0.1r)'‘  1 

L  fc'  J 

fc=2,“ 

=  (1.5  X  10-2,8.9  X  10-4,3.9  x  lO-^,. 

As  can  be  seen  the  error  is  small  and  has  its  largest  values  at  the  edge  of  the  aperture. 

In  order  to  get  a  general  estimate  of  the  error,  assume  Vj,  k  >2,  and  we  obtain 

£:(r)  <  1.5  X  10-2  ,/,(2)(aj.)  (34) 

For  (j)^'^\ar)  «  (j)^^\ar),  the  error  is  <  1%,  hence  small  but  still  measurable.  The  table  below  demonstrates  how 
this  technique  works  in  the  presence  of  a  sinusiodal  distrubance  of  amplitude  tt.  The  leftmost  column  is  the 
radius  of  the  pupil,  the  second  column  the  amplitude  of  the  phase  disturbance,  the  third  column  the  magnified 
pupil  coordinates,  the  forth  column  the  amplitude  of  the  magnified  phase  (the  beacon  phase),  the  fifth  column 
the  correction  term  as  given  in  Equation  32,  the  sixth  column  the  corrected  beacon  phase,  and  the  last  column 
the  error.  As  can  be  seen,  for  this  example,  the  technique  works  as  calculated. 


Application  of  the  technique  to  a  test  phase  of  7rcos(r).  Magnification  =  1.1  | 

radial 

distance  (m) 

science 
beam’s  phase 

Correction 
—  A(j()(r) 

magnified 
rad.  dist.  (m) 

beacon 
beam’s  phase 

corrected 
beacon  phase 

error 

0 

3.14 

0 

0 

3.14 

3.14 

o.ooE-eooo 

0.25 

3.04 

-0.02 

0.28 

3.02 

3.04 

9.49E-004 

0.5 

2.76 

-0.08 

0.55 

2.68 

2.75 

3.41E-003 

1 

1.7 

-0.26 

1.1 

1.43 

1.69 

8.04E-003 

1.75 

-0.56 

-0.54 

1.93 

-1.09 

-0.55 

1.13E-002 

Wave  Optics  Simulation  of  the  Correction  In  order  to  test  the  correction  in  a  more  realistic  environment, 
a  wave  optics  simulation  was  run.  The  simulation  used  WaveProp^^  to  create  recreate  the  DNiP  sodium  guidestar 
geometry.  Typical  results  are  plotted  in  Figures  6  and  7. 

Figure  6  displays  phase  for  a  single  frame  of  a  single  realization.  The  upper  left  plot  is  the  phase  of  the  science 
beam  corrupted  by  the  atmosphere.  The  beason’s  measured  phase  is  displayed  in  the  upper  right  plot;  this  is 
the  phase  that  is  used  as  the  estimate  of  the  science  beam  phase,  i.e.  it  is  the  phase  placed  on  the  DM  which 
corrects  the  science  beam.  The  middle  left  plot  is  the  correction  calculated  from  the  beacon  phase  to  be  added 
to  the  beacon  phase.  The  middle  left  plot  shows  the  beacon  phase  corrected  via  Equation  32.  This  corrected 
beacon  phase  is  placed  on  the  DM.  The  bottom  plot  displays  a  radial  average  of  the  residual  error  before  and 
after  correction  (recall  in  simulation  we  have  access  to  the  science  beam  phase  so  it  is  used  as  truth  phase).  The 
blue  curve  is  the  uncorrected  phase  which  shows  high  residual  phase  error,  and  the  red  curve  is  the  corrected 
phase  showing  greatly  reduced  residual  phase  error.  Note  that  the  degradation  is  so  small  to  be  unobservable  to 
the  eye  in  the  2-D  plots. 

Figure  7  displays  phase  for  a  different  realization.  The  top  row  is  the  phase  of  the  science  beam.  The  middle 
left  plot  is  the  residual  phase  error  if  this  phase  is  corrected  using  beacon  phase  measured  in  the  DNiP  plane. 
The  middle  right  plot  is  the  residual  phase  error  if  this  phase  is  corrected  using  the  corrected  beacon  phase. 
The  bottom  row  displays  a  radial  average  of  the  residual  error  before  and  after  correction.  The  blue  curve  is 
the  uncorrected  phase  which  shows  high  residual  phase  error,  and  the  red  curve  is  the  corrected  phase  showing 
greatly  reduced  residual  phase  error. 


7.  DISCUSSION 


There  remains  much  work  to  do.  For  instance,  we  don’t  discuss  how  to  estimtate  the  derivative  in  the  presence  of 
beacon  noise  or  how  often  should  the  gradient  reference  be  updated.  Also,  in  the  dynamic  case  one  must  use  the 
DM  commands  to  estimate  the  beacon  phase,  i.e.  use  (('experimentai(ar)  =  0dm  +  ^residual  to  calcualte  0*-^^(ar), 
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and  we  don’t  discuss  implications  of  DM  hysteresis  in  this  estimation.  Also,  we  don’t  discuss  loop  stability  in 
the  closed  loop  case  (recall  by  adding  a  time  varying  reference  to  the  gradient  map  we  are  effectively  adding  a 
additional  small  feedback  loop  into  the  closed  loop  AO  system) . 

However,  we  did  show  that  the  degradation  in  performance  of  the  proposed  SOR  high  efficiency  AO  design 
is  due  to  the  DNiP  geometry  and  not  caused  by  focus  anisoplanatism,  tq,  9q,  or  We  also  demonstrated  that 
the  culprit  for  the  DNiP  degradation  is  magnification  and  given  an  expression  to  both  describe  the  degradation 
and  a  theoretical  means  to  mitigate  the  degradation. 
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APPENDIX  A.  RELATION  OF  THE  SCIENCE  BEAM  TO  THE  BEACON  BEAM 

First  note,  the  residual  focus  error  at  the  pupil  plane  is  given  by  e  ^  ^  From  the  lens  equation, 

A  -I-  L  =  -1—^  Since  /  =  Zsd, B,nd  Zim  is  the  effective  virtual  image  location  of  the  expanding  beam  on  the  right 
side  of  the  lens.  So  the  residual  focus  error  is  ^  ~  When  the  propagation  is  extended  to  the 

DM  plane,  then,  the  new  residual  focus  error  in  the  DM/WFS  plane  is  where  R'  =  R-\-  Zdm,  so 

j_  _  1 
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Suppose  there  is  some  field  Umeas{x,  y)  at  the  DM/WFS  plane,  measured  by  the  WFS  focused  at  the  beacon. 
This  is  the  field  measured  by  the  WFS,  and  there  should  be  no  residual  focus  error  on  it.  We  are  going  to  use 
this  field  to  re-create  the  field  from  the  beacon  seen  by  the  DM.  When  the  residual  focus  error  at  the  DM/WFS 
plane  is  added,  the  field  becomes: 


U{x,y)  =  17meas(a:,  j/)e’2 


(35) 


Back  propagate  this  to  the  pupil  plane.  (For  back  propagation,  k  =  -k). 


^—ikzdr. 


Ubacki^i  — 


iXz, 


dm  J  JWFS 


dxdyUmeas{x,y)e  ^ 


fc /-...s  1  ,,2n  1  i 


'-§(x^+y 


I  ^  g  Zzd„ 


(36) 


Now,  remove  the  residual  focus  error  at  the  pupil  plane  (effectively,  send  it  through  the  lens  with  /  =  Zsci)- 
Then, 


=  ^Uback{^,y) 

Finally,  propagate  back  to  the  WFS/DM  plane  and  compare  to  the  original  Umeas(x,y)- 

aikZdr, 


U{x',y')  = 


'^XZdm 


d^drje^ 


-iii-x'f+iv-y') 


(37) 


(38) 


pup 


Now,  simplify  this  equation  to  obtain  a  relationship  between  U{x',y')  and  Umeas{x,y).  The  equation  for 
U(x’,y’)  with  all  terms  included  is  shown  below 


U{x',y')=- 


oikZdr, 


iXZflm 


-ikZd-n 


J  JWFS 


For  now,  we  can  work  with  the  ID  case  U{x').  (The  2D  case  has  y,  rj  and  y’  added  in  the  exponents  to  y, 
and  x’,  respectively,  so  it  is  trivial  to  extend  to  the  2D  case  after  simplifying  the  equation). 


U{x')  = 


o'^kZd-n 


-ik  c2  g  '^kZdr, 


dxUmeas{x)e'^'^  R'e^^dr 


(39) 


'WFS 


To  simplify  this,  put  everything  inside  all  of  the  integrals  and  expand  the  quadratics  in  the  exponents.  Group 
terms  with  like  x,  y,  x’,  y’,  rj,  and  take  x’  and  y’  terms  out  of  the  integrals.  Then  switch  the  order  of  integration, 
which  is  possible  because  the  fields  are  smooth  and  finite. 


U{x')  = 


1 


/  d^  dxUmea8{x)e’‘d^^ 

I  pup  JWFS 


,  1  Zk  /2 


dxUmeas{x)e’^'‘  ^ ^  /  dS,  e 


'WFS 


(41) 


/  pup 


Look  at  the  inner  integral,  over  the  pupil.  By  completing  the  square,  ax^  +  bx  +  c  =  {^/ax  +  \fc)^ ■  Then 
b  =  2^/ac.  So 


1  ,  Ax- x')  1  (a; -a;') 


=  (V5=f 
^dm  V  ^^sci 


^dm 


-^sci  \2 

o  J 


{x  -  x'Y  Zs 


(42) 


‘^dm 
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^ \f^^-  Then  the  integral  is 


^/2^ 


^  {x-x')^  Z^ci 


du 


—  iku^ 


>  pup 


(43) 


This  integral  is  known  to  be  — -L-U — i) — Vku)  ^  recall  that  (—1)  = 


i2y/k 

■/2i-e  ..  [ - - 


u=  — oo 


(44) 


The  error  function  gives  1  and  oo  and  -1  at  —  oo,  so 


's/'^^sci' 


ik- 


2Vk  2Vk  V  k 

Throw  this  back  into  the  equation  for  the  field  at  {x',y'),  and  you  have: 


TTf  \/2'nZsci  **=  f  j  TT  f  \  — 

[/(x)  =  ^=— - —e^-drr,  /  dxUmeas[x)e  2-d™^e 

Vk{\Zdm)  JWFS 


(45) 


(46) 


Expand  the  quadratics  again  and  simplify: 


,  V^TTZsci  f 

Vk(XZdm)  JWFS 


,  ,  -7-9.2C  1  1  ■>  ik  “V'-  (x^  —2x' x+x'‘^) 

dxU^eas{x)e  -^dJe  ^‘d^' 


TT/  l\  '\/2FZsci  *^2;  (2*1  +2/2^’  1  /"  77-7-  /  \  ^  H'  2zl  +2zT*  ^  a;(  x  ) 

[/(x)  =  - -e  ""rf™  /  dxUTneas{x)e  ” 

Wk[\Zdm)  JWFS 


(47) 

(48) 


Let  — \-  Complete  the  square  again,  with  a  =  x/a^  h  =  c  =  ■  Then 


TT/  /\  ‘sJ'^'^^SCi  (2^^™  + 

LlX  )  =  — ;= - e 

^{\ZdmY 


dm  22“  4 


dm  dm  I  dx  Urneas  (^X^  C 
JWFS 


iKVax+  f 


(49) 


This  can’t  be  analyzed  further  without  Umeasix).  Let  ^  =  55^ — h 


.  Then  the  2D  version  looks 


like  this: 


U{x',y')=  f  f  dxdyUmeas{x,y)e 

WkyXZdm)  J  JWFS 


ik{{y/ocx-\-  2  +(v^y+ 


y  ^sci  \2^ 


(50) 


Any  possibile  mitigation  of  the  degradation  only  considers  phase,  so 

Arg{U{x' ,y'))  =  Arg  [-^{x'"^  +  y'"^)  [  [  dxdy  Umeas{x,y)e 
\z  J  JwFS 


(51) 
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Residual  Phase  -  radial  average 


Figure  6.  Top  row:  science  beam  phase  (left)  and  beacon  beam  phase  (right).  Middle  row:  correction  phase  (left)  and 
corrected  beacon  phase  (right).  Bottom  row:  radial  average  of  residual  phase  error,  blue:  uncorrected  phase  showing  high 
residual  phase  error,  red:  corrected  phase  showing  greatly  reduced  residual  error.  As  can  be  seen  the  degradation  and 
the  correction  terms  are  small  but  measureable. 
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Figure  7.  Top  row:  science  beam  phase.  Middle  row:  residual  phase  error  uncorrected  (left)  and  corrected  residual  phaser 
error  (right).  Bottom  row:  radial  average  of  residual  phase  error,  blue:  uncorrected  phase  showing  high  residnal  phase 
error,  red:  corrected  phase  showing  greatly  reduced  residual  error.  As  can  be  seen  the  degradation  and  the  correction 
terms  are  small  but  measureable. 
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